SURVEY ON MODERN RADAR 
SIGNAL PROCESSING 



Joao Paulo Goncalves Barcia 



OUOLEV KNOX UBBARX 
xAVAL POSTQB*i>UATE SCHOOL 
MONTEREY. CAUFORNIA 93M0 





>^22a*v>, ^ 



i ii r>.i;sf9£3 






mraeooBkBttxvm 
! ^ 



««*wcsra 



SURVEY ON MODERN RADAR SIGNAL PROCESSING 



by 



Joao Paulo Goncalves Barcia 



December 1975 



Thesis Advisor: 



John Bouldry 



Approved for public release; distribution unlimited. 



T171667 



Unclassified 



SECURITY CL ASSlFICATiOR OF THIS PAGE (V'^han Dmtm Entorod) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPOFiT NUMBER 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


4. TITLE 

Survey on Modern Radar Signal 
Process ing 


5. TYPE OF REPORT » PERIOD COVERED 

Master's Thesis; 
December 1975 


6. PERFORMING ORG. REPORT NUMBER 


7. authorc*; 

Joao Paulo Goncalves Barcia 


8. CONTRACT OR gran T number^'*; 


S. PERFORMING ORGANIZATION NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


10. PROGRAM element. PROJECT, TASK 
AREA ft WORK UNIT NUMBERS 


1 1. CONTROLLING OFFICE NAME AND ADDRESS 

Naval Postgraduate School 
Monterey, California 93940 


12. REPORT DATE 

December 1975 


13. number of pages 


t4, monitoring agency NAME & ADDRESSf/f /rem Controlling Oflico) 

Naval Postgraduate School 
Monterey, California 93940 


15. security class, (of thia reporO 

Unclassified 


15*. DECLASSI FI cation/ DOWN GRADING 

schedule 


16. DISTRIBUTION ST ATEMEN T fo/ fh/# K^porO 

Approved for public release; distribution unlimited. 


17. DISTRIBUTION STATEMENT (oi tha ebetfct cn/crad In Block 20, If dliforent from Report) 


18. SUPPLEMENTARY NOTES 


IS. KEY WORDS (Continue on reveree aide If nacecamry and Idantify by block number) 

Radar 

Digital 

Signal processing 



20. ABSTRACT (ContInuB on rBVBfao aidm H nmcmaoAry mnd identify by block numbat) 



The purpose of this thesis is to investigate the state of 
the art of radar signal design as well as radar signal processors 
and determine the actual trends in modern radar design. The use 
of a digital general purpose radar signal processor is discussed. 
The concepts of ambiguity and autocorrelation function are 
investigated in regard to radar resolution capabilities. The 
concept and analytical development of the DFT/FFT is presented. 



DD 1473 

(Page 1) 



EDITION OF i MOV «5 IS OBSOLETE 
S/N 0102* 014- 6601 I 






SECURITY CLASSIFICATION OF THIS PACE Data Bntarad) 



1 



Unclassified 

S‘lrCOKl'’^V CLASSIFICATION OF THIS PAGEfH^en Dft« Ent0f<ij‘ 



Quantization noise in a digital MTI processor and its effects 
in the improvement factor are analyzed. Optimization techniques 
for the response curve of digital MTI processors using staggered 
PRF are investigated. The SAR concept and analysis as well as 
techniques to obtain low correlator rates in SAR digital 
processors are presented. 



DD Form 1473 
, 1 Jan 73 

S/N 0102-014-6601 



Unci a ss i f.ie d 

security CLASSiriCATION OF this PAGEf»7i»n Enlrrtd) 



2 



Survey on Modern Radar Signal Processing 



by 



Joao Paulo Goncalves ^arcia 
First Lieutenant, Portugtiese Navy 
B.S., Naval Postgraduate School, 1974 



Submitted in partial fulfillment of the 
requirements for the degree of 



MASTER OF SCIENCE IN ELECTRICAL ENCINEERINC 



from the 



NA^/'AT DnCTT'n A TMT A 



I Zies/s 

6s./ n<^ 



PUDLEV ICJOX LISRAR* 

NAVAL postgraduate SCHOOL 
MONTEREY. CAUFORNU 83840 



ABSTRACT 



The purpose of this thesis is to investigate the state 
of the art of radar signal design as well as radar signal 
processors and determine the actual trends in modern radar 
design. The use of a digital general purpose radar signal 
processor is discussed. The concepts of ambiguity and auto- 
correlation function are investigated in regard to radar 
resolution capabilities. The concept and analytical devel- 
opment of the DFT/FFT are presented. Quantization noise in 
a digital MTI processor and its effects in the improvement 
factor are analyzed. Optimization techniques for the response 
curve of digital MTI processors using staggered PRF are in- 
vestigated. The SAR concept and analysis as well as techniques 
to obtain low correlator rates in the SAR digital processors 
are presented. 
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I. INTRODUCTION 



The processing of signals has been always a necessity 
whenever information was to be transmitted through a channel 
between a source and a destination. It was not until 1948 
that Shannon presented the basic theory of information re- 
lating the source entropy with channel capacity and the 
probability of error. Techniques of signal processing were 
already in use, but the results given by the information 
theory clarified the limits that could be achieved. 

Early analog processors with enough bandwidth and reason- 
able signal to noise ratio could process data at extremely 
high data rates, but the components' low stability and high 
cost, and sometimes low versatility, were responsible for a 
gradual substitution of analog by digital processors when 
the price of digital logic went down, and the digital hard- 
ware technology suffered a big jump in the last ten years. 

In fact, in today's radar signal processing the trend to 
a substitution of analog by digital processors is inevitable. 
The theory and boundaries of digital signal processing are 
continually being improved as well as hardware components. 

Speed is almost no longer a problem and versatility is immense. 

Chapter II is an introduction to radar history. In Chapter 
III signal processing is related with the radar problem through 
the radar range equation and a general digital processor for 
a radar system is presented. In Chapter IV an analysis is 
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made of the present signal processing techniques in radar 
with special emphasis into digital approaches. In Chapter V 
an investigation is made into the recent developments and 
problems of tv\fo major areas of radar digital signal process- 
ing. In Chapter V an overall view of the achievements and 
trends of radar signal processing is explored. 
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II. RADAR. HISTORY AND APPLICATIONS 



The electromagnetism and electromagnetic wave propagation 
theories of Maxwell, later (1886) experimentally verified by 
Hertz, contain all the necessary background to understand the 
principles of radar. The first experiments on detection with 
radio waves started in 1903. Due to inadequate technology 
the results were very poor. The obtained ranges were less 
than those achieved with optical systems and this was enough 
to discourage any interest in pursuing the experiments. It 
was only in 1922 at the Naval Research Laboratory that a 
wooden ship was detected using a CW radar with the transmitter 
and the receiver as separate units. From then on an increas- 
ing interest in radar technology became evident. In 1930 
using a bistatic radar, the first aircraft was detected, and 
by 1932 the detection ranges were already in the 50-mile 
region. By that time the pulse techniques were not yet 
explored, so the information was in the presence or absence 
of a target and not in range information. But in 1935 the 
British successfully used pulse techniques to measure distance 
to the targets. In 1936 the detection range was already in 
the 90-miles region. The basic principles were understood. 

Its next steps were in technology improvements in order to 
get stable elements and more power output to increase the 
range coverage. Higher frequencies were used in order to 
work with smaller elements and high gain antennas. 
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The discovery of the magnetron was an important step since 
powers went up by a factor of 100, and wavelengths of 10 cm 
were obtained. But the greatest incentive for further devel- 
opments in the radar field came from the military necessities 
of World War II. After the War a stagnation of about five 
years slowed down the rhythm of improvements in the field. 

But in the fifties the introduction of the high power klystron 
increased not only the power output but the frequency stability 
necessary for coherent detection as well as coherent MTI ; also 
very low noise receivers were implemented. 

From then on, with more and more perfect technology, more 
complex systems appeared. The introduction of the computer 
as a storage and control element was of extreme importance. 

All ballistic systems and very large coverage systems suffered 
a big impulse. Digital systems improved the implementation 
of stereable array antennas, synthetic aperture radar, track- 
ing systems, as well as all digital processing techniques. 

With the development of radar, the statistical nature of 
the detection problem gave rise to an increasing interest in 
the study of the statistical properties of clutter, radar 
cross section, rain, etc. So, due to a more precise knowledge 
(modeling) of the systems environment, the radar systems 
changed depending on the type of application. 

Thus, a design to optimize a given application seldom 
can be used in others. Today the radar design field is wider. 
Applications vary from military to civilian to scientific. 

With improvements in technology and more accurate modeling, 
research in the field is far from becoming saturated. 
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III. SIGNAL PROCESSING IN RADAR 



A. INTRODUCTION TO THE RADAR RANGE EQUATION 
Only the analysis of the equation 

R = £(nj^>n2> ••• (3-1) 

where R is the distance from the radar antenna to the target, 
and (nj^,n 2 > ... n^) are the radar and environmental parameters, 
can give the necessary knowledge to maximize R. If is the 
peak pulse power radiated and G the gain of the antenna, the 
transmitted power density is 

Ai 

4 TT R^ 



and the received power density, after reflection in a target 
of cross section will be 



Pt" ^t 

C4itR2) 2 



(3-2) 



If (3-2) is multiplied by the antenna aperture A^ , the 
received power will be 



P 



r 



PtG cjj 



(4itRA 



2 -i 2 



A 

e 



(3-3) 



If is the minimum received signal, 

distance is 



R“ 

max 






(4it) ^S 



min 



then the maximum 



(3-4) 
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(3-5) 



Since by definition of noise figure (Fj^) of a receiver 



S. = N. F ( -j^) 
1 1 n '' Nq 



i = input 
o = output 



where 



N. = kTB = input noise 
1 n ^ 



(3-6) 



Then, combining (3-4) with (3-5) and (3-6) 



R" 

max 






(4,r)^kTB F„(S,/N.) . 

n ri'- O' O'* mi 



in 



Since in a pulse radar the average power is 



P = T f 
av t r 



T = pulse width 
f^;. = pulse rejection frequency 



and the receiver bandwidth is approximately 1 /t 



R“ 

max 






(4Tr)'kT£ F (So/N„) . 

^ r n^ ° o-^min 



(3-7) 



where (S /N ) . is the minimum signal to noise ratio at the 

output of the linear section of the receiver, necessary for 
detection. The first important characteristic in equation 
(3-7) is the fact that some of the variables can be considered 
deterministic but others must be analyzed as random variables. 
The random variables are and S/N. Statistical discriptors 
for are already extensively studied [1-2] not only in the 
general case but also for some specific cases, as for different 
types of clutter. Since the noise is normally treated as 
Guassian, all necessary statistical discriptors are available. 
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The statistical discriptors together with the decision 
criteria will determine the probability of detection and 
probability of false alarm. 

B. RELATION BETWEEN THE RADAR RANGE EQUATION 

AND SIGNAL PROCESSING 

After the War (1945) most of the effort in increasing 
range capabilities (eq. #3-7) was focused on increasing 
average povv^er (P ) , antenna gains (G) , and mixers and 
receivers with lower noise figures (F^^) • With the magnetron 
and the klystron, the peak powers reached the megawatts range. 
Soon it was verified that since the technology was already in 
the limits of transmitted peak power, any increase in peak 
power was questionable from a technical and a financial aspect. 
Also, since 

R = kP , k = constant 

In R = In k + i In P 

dR = 1 ^ 

R 4 P 

That means, that the percent gain in range is only 1/4 of the 
percent gain in power. The pulse width could not be increased 
much more because of the required range resolution. This led 
to the development of what today is called signal processing 
techniques. Waveforms were studied in order to increase the 
average power but maintaining range resolution. The basic 
teclmiques will be presented in Chapter IV. 

The brute force approach to increase the range of detec- 
tion. (3-7) was used until 1955. The technology then shifted 
to an improvement of the signal to noise ratio. If the new 
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signal to noise ratio is 

(S/N) = D(S^/N^), D > 1. 

Since D > 1, there is an increase in the maximu... rang/ of 
detection (3-7) . 

C. USE OF COMPUTERS AS SIGNAL PROCESSORS IN RADAR 

As range and aximuth resolution became finer, as multi- 
target processing due to digital techniques became more 
effective, the conventional final processor in the radar 
system, the man, gradually was substituted by the computer. 

The total cimount of information and the necessary speed for 
processing, tax the capabilities of man. So, digital tech- 
niques are gradually being substituted for not only those 
processors that used analog devices, but also where human 
operators were used. 

Figure 1 gives an overview of the techniques presently 
used in digital signal processing. As in the analog theory 
linear approximation is used. So, using the theory of linear 
discrete time invariant systems there are two common approaches 
digital filtering and spectral analysis. Using the theory of 
digital filtering, algorithms can be implemented that correspond 
to finite impulse response filters (FIR) or infinite impulse 
response filters (IIF). Using spectral analysis, two methods 
are possible: implementation of Fast Fourier Transform (FFT) 
algorithms or use of statistical spectrum analysis. 

Both digital filters and spectral analysis are affected 
by quantization noise, which must be taken into account. 
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As seen in Fig. 1, the radar appears as one of the prime 
applications for digital systems. 

It was in the 1960 's that the first ]ow data rate digital 
processors were implemented in radar systems, basically to 
perform tracking and weapon order computations. The next 
step was the implementation of the MTI processor. A flexible 
programmable digital processor was not achieved until very 
recently. The low cost of semi-conductor memory and the 
advanced technology in microprogramming, witli inherent increase 
in speed, were the main reasons for this achievement. Now a 
typical general purpose signal processor for radar will be 
described. Figure 2 shows a block diagram of a digital signal 
processing system for a radar. The main processing unit is 
the arithmetic pipeline, which is a general purpose signal 
processor (GPSP) . This consists of five sections selected 
in a way that the most frequently necessary operations can 
be performed; Fast Fourier Transform (FFT) , recursive and 
nonrecursive digital filtering such as MTI, threshold gener- 
ation, peak detection noncoherent video integration and range 
and angle estimation. The first block in the pipeline, the 
matrix switch, selects data from data memories, also from the 
receiver through a track buffer, or from anoth er GPSP via a 
wrap-around connection. The second block, data scaling, 
connects the data for floating point arithmetic. The complex 
multiplication block does the equivalent for four real multi- 
plications and two real additions, which corresponds to a 
complex multiplication. Since some operations, like integration. 
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ARITHMETIC PIPELINE 
(GPSP) 




The re-normalization block is sometimes needed. Connected 
to the arithmetic pipeline is a memory section and a control 
section. The output of the arithmetic pipeline is relayed 
to an input/output control unit (IOC) which is a buffer for 
the radar data as well as a macro control for the processing 
functions of the GPSP. The advantages of this processing system 
are that: once the data is started the process is continuous, 
and does not follow the normal sequence of fetch, operate, and 
store, characteristic of general purpose computers; the pro- 
cessor has a horizontal structure as opposed to the normal 
vertical structure of computers; this means that an instruc- 
tion, after executioii, causes all operations specified to be 
executed in a clock cycle; there is a hardware separation of 
instructions and data in the memory, which simplifies the job 
of the programmer; the control memory is made of RAM's which 
implies a much greater flexibility in subroutine changes; and 
the execution of instructions is w'ithin the same time frame 
that corresponds to a range or doppler cell. 

Four memories interact with the arithmetic pipeline. Two 
of them store and recirculate data, the third stores weighting 
coefficients for such filters as MTl, FFT or pulse compression, 
the fourth stores micro opcodes. Data sequences and arithmetic 
operations are both controlled through software, by two control 
units. The program that resides in the sequence control memory 
(SEQ- CONTMEM) , informs the machine of range and doppler dimen- 
sions of the problem. The program in the Arithmetic Control 
Memory, ACM, informs the machine of the particular algorithm 
to be used. 
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IV. VARIOUS ASPECTS OF THE THEORY 
OF RADAR SIGNAL PROCESSING 



A. SIGNAL PROCESSING 

The desire to transmit information is closely related to 
the processing of signals. Since normally the information to 
be transmitted cannot flow through the available channels, it 
is necessary that it be processed to use the chosen channel. 
Looking at Fig. 3, the inverse process has to be accomplished 
in order to get the information in an understandable form at 
the receiving end. The process is not so simple, since in 
general the processors themselves and the channel introduce 
noise that may or may not eliminate or change the understanding 
of the message. So there is a need to interpret the results 
with some decision criteria. Here decision theory plays an 
important role. Even for the same type of application, for 
instance, radar, decision criteria are not the same. The 
criteria may be of constant false alarm rate, or may be of 
fixed probability of detection. Basically every application 
or system must have its decision criteria. 

There are two main approaches to the processing of signals, 
the time domain approacli and the frequency domain approach. 
Besides the fact that they use completely different hardware, 
in the pure mathematical way, they are duals of each other 
through the Fourier transform pair. So the only reason for 
one type of approach is due to more perfect reliable type of 
hardware . 
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With the advance of computer techniques , both approaches 
are split between digital and analog signal processing. 

It was only after World War II that electronic engineers 
became interested in the applicability of digital hardware 
techniques in signal processing areas. But it was not before 
the late sixties that digital techniques suffered a big jump 
replacing some of the analog techniques. The foundations of 
digital signal processing can be related to tlie Laplace's 
Z- transform theory; but only in the mid sixties was the theory 
presented in a formal way. Various papers appeared at that 
time related to this subject, but only in 1969 was the first 
attempt made for a comprehensive theory of digital signal 
processing [21]. Most recently (1975), two very comprehensive 
books by Oppenheim [3] and Rabiner [4) can be considered as 
giving an excellent treatment of the subject, not only through 
its mathematical structure but with very good applications. 
Today's trend is definitely towards the substitution of analog 
processing by digital processing. The speed achieved in 
today's digital processors, the very fast algorithms used to 
implement the FFT, the gradual substitution of infinite impulse 
response metliods to finite impulse methods due to a better 
knowledge of theory and higher efficiency of calculations are 
the major factors responsible for the increasing shift to 
digital techniques. There are yet some type of applications 
where the speeds are so high, or the digital hardvjare is so 
complicated that analog techniques are still used. 
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B. AMBIGUITY FUNCTION AND AUTOCORRELATION 



1 . Concept 

Considering only non- accelerating targets, the two 
kinds of imf ormation of interest are the position of the 
target and its velocity. Given two targets, the capacity 
to differentiate between their positions and velocities is 
vital for most radar applications. 

The capacity to differentiate in range is called 
range resolution, and in velocity, velocity resolution. 

As far as range resolution is concerned, it is 
obvious the shorter the pulse the better is the differentia- 
tion between the two targets, that is, the higher is the 
range resolution. With the measurement of velocity, since 
it is directly related to the phase difference between sig- 
nals, the longer the pulse the greater the number of cycles 
that can be compared and consequently the better the doppler 
and the velocity resolution. So, at first glance it looks 
like to optimize one of the resolutions the other will have 
to be sacrificed. The other solution is to try to get an 
optimum solution for both cases at the same time. There is 
a quantitative way to express these conditions in a precise 
way. That was the reason for the appearance of the mathe- 
matical concept of ambiguity function. 

2 . Range Ambiguity Function 

Let the transmitted signal be represented by S(t) = 
Re[i|;(t3], i[)(t) = u(t) exp. jw^t. Then, with no doppler 
present the received signals from two different targets are 
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H'(t) and 

T= equivalent time difference between the two targets. 

If it is chosen, as a measure of the resolution, the mean square 
of the difference between the two received signals, then 

+ 00 

e^(t) = - i/j(t-T)|Mt 

= rdt + _/Ji|;(t - t) I ^dt 

+ 0O 

■ (t)ij;(t - T)]dt 

-Hoo 

= 2 /_^ |u(t) I ^dt - 2Re[exp -j(jj^r[/u*(t)u(t - x)dt]] 



Since the first part of the result 

total energy of the signal and has 

second term is going to make e^t) 

+ 00 

= loo u*(t)u(t - x)dt. 



is proportional to the 

a constant value, only 

vary. Therefore 

* = complex 
conjugate 



the 



is defined as the range ambiguity function; only c(x) varies 
the value of e^(t). As can be seen, there is a perfect 
identity between c(x) and the autocorrelation of the signal. 
The ideal case would be to have the ambiguity function with 
a spike at the origin and zero anywhere else. That would mean 
that the only situation where it was impossible to differen- 
tiate the two signals was when they were on top of each other. 

As an example. Fig. 4(a) shows the autocorrelation 
function of a pulse of duration T, Fig. 4(b) its uncertainty 
function ]c(x) |^. In part (c) there is a rectangle with the 
same height of c^(o) and with the same area under |c(x)|^. 
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that is, A = B. As it can be seen, the base of that rectangle 
is a measure of the spread of the curve and is called the 
delay resolution constant; so 



At 



_ /^“|c(t) |Mt 
“ - 00 “ 

C2 (o) 



sometimes is more practical to define At as a function of the 

+ 00 

bandwidth. Let u(to) = F[u(t)] = u(t) exp - jwt. 

since by Parsevall's Theorem F [c(t)] = |u(m)|^ 

2tt I u (w) I "* dm 
At = 

[/^"“luCm) Pd'03]^ 

- 00 



if the effective bandwidth is defined as 



[/.';;iu(a3) |2dco]2 

4tt/^'”|u(w) I '*dw 



At = 



2 W 



and the range resolution becomes 



A R = 



cAt 

2 



All these definitions make sense since as pointed out 
before, with a short pulse there is a higher ran^ resolution, 
but to have a short pulse there is need of a wider bandwidth, 
so there is equivalence in the two expressions for At, 

3 . The Velocity Ambiguity Function 

Using the same type of mathematical logic and defining 

i|j(f) = F[(})(t)] 
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a complex correlation function is determined [5] to be 



K(f^) = /u*(2TTf)u(27rf - 2-rrf^df 
= /|u(r)|^ exp j2ir£^tdt 

and the doppler resolution constant [5] 

Af = /|K(fd) IMfd ^ /|u(t) rdt , J_ 
kHo) [/|u(t) 1 2dt]2 Te 



where Tg is the effective duration. 



Since f^j = 



dfd = 



2vf o 
c 

2dv fo 



which infers Afj^ = - 



2 fo 



A V 



the velocity resolution constant is 
c 



Av = 



2 fo Te 



The graphical interpretation for Te is similar to the 
graphical interpretation given for Ax in Fig. 4 but instead 
of using the function |c(x)|^, |K(fj)|^ is used. But as seen 

earlier, the optimization of one resolution implies the mini- 
mization of the other. This implies the need to study a two 
dimensional correlation function in order to see the mutual 
effects between the two resolution parameters. 

So if the transmitted signal is 

i|;(t) = u(t) exp j2Trfot 

the received signal with doppler and delay will be 
i{;(t-T) = u(t-x) exp j 2TT(fg - fd) (t-T) 
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The mean square of the difference 



will yield a two dimensional correlation function of the form 

X(r,f(j) = /u(t) u*(t-r) exp j27rfjtdt 
Easily it can be seen that 
X(T,0) = C(T) 

X(0,fd) = k(f^) 

If the volume under |X(T,f(j)|^ is determined and 
divided by |X(0,0)|^, an equivalent resolution parameter 
called the effective area of ambiguity is obtained. 

)= /-^l XCt , f^) I ^dxdfd 
|X(0,0)|2 

It is easy to prove [6] that 

//|X (T,fd) l^dxdfd = |X(0,0) 

A('T>fd) = 1 

So if an effort is made to improve one dimension, the other 
will never get better. But given one acceptable resolution 
in one dimension, there are ways to find the optimal solution 
for the other dimension if X(r,f^) is known. Fig. 5 is an 
ambiguity surface for a monochromatic pulse with a smooth 
envelope. It would be a reasonable temptation to try to 
design the waveform from a desired ambiguity function. That 
turns out to be impractical. First, a criteria to design the 
ambiguity function would be difficult to build. Secondly, and 
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Fig. 5 
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most importantly, when designing a radar system normally there 
are other more important factors, from economic to space and 
weight constraints as well as technical, that dictate the 
boundaries for the waveform that we have to pick. Also, there 
is an important relation between the chosen waveform and the 
type of clutter model , and the approximation of the clutter 
model. Normally some flexibility in the digital processing 
is desirable. 

C. PULSE COMPRESSION 
1 . Concept 

Soon it was found that trying to increase the maximum 
range by increasing power, keeping the same pulse width, v/as 
the difficult way. Isolation problems were difficult to 
handle as well as reliable components. The only way to in- 
crease the energy of the signal was by increasing the pulse 
width. But an increase in pulse width would reduce the range 
resolution. Since the real problem was to increase the time 
of transmission but at the same time do not decrease the 
bandwidth (effective) of the waveform some different types of 
waveforms were studied, using the ambiguity function concept. 
The processing of the signal should be such that the trans- 
mitted waveform with high time of transmission and small 
effective bandwidth (Wg) should be converted to a high energy 
pulse with pulse width approximately 1/Wg. This is called 
pulse compression. There are two basic ways of implementing 
pulse compression techniques. The active way, which is a 
time domain approach which basically uses correlators in the 
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detection and active elements to produce the waveform in the 
transmission. The passive way, which is a frequency domain 
approach, uses passive filters to generate the waveform and 
matched filters in the receiver. Each of the two has some 
variations, and it is even possible to build a system as a 
combination of the t\\’o processes. Table I gives a relative 
performance of various types of pulse compression techniques 
2 . Linear FM (Chirp) Pulses Compression 

First, the passive generation of linear FM signals 
will be analyzed. An IF pulse generator feeds a dispersive 
delay filter with a frequency versus time characteristic as 
in Fig. 6(c). The signal is then up converted and only one 
of the bands is transmitted. At the receiver the signal is 
down converted to IF and then passed through a dispersive 
filter with the opposite slope (Fig. 6 c) . Using at recep- 
tion a mixer sideband inverter, the same type of dispersive 
filter can be used because as it can be seen, £2 and f^ 

(Fig. 6b) have opposite frequency versus time character- 
istics. In fact, with this technique the same dispersive 
filter for both transmission and reception can be used. 

The analysis approach to this type of processing is 
the following: 

Let the received signal be of the form 
f(t) = exp j[(oJo + W(i)t + ^yt^], |t| < T/2 

= 0, |t| > T/2 

where y = = frequency slope of the dispersive filter 
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T = duration of the transmit pulse envelope 
o)(i = doppler shift. 

The pulse spectrum is 

+ 0O 

F(to) = exp -jmt dt 

+ 00 

= -Coo j[(fjJo + 03d - m)t + 

The filter transfer function is 
H(co) = exp j [ (wo - w) V2y] 

which implies G(m) = H(m)F(o3) = exp j [ (wq - a3)^/2y] 

T/2 ^ 

J exp j [ (wo + cjd ■ ^ T F 

-T/2 ^ 

which implies the time function at the output of the filter 
is 



g(t) = F'^[G(m)] = ^ / ^G(to) exp jmtd(^ 



+ 00 



which implies g(t) = f ^ exp j (wq ' m)^/2y 



T/2 ^ 

f exp j [ (wq + md - w)t + yyi^] 
-T/2 



dr exp jcjt dco 

Inverting the order of integration: 

T/2 +00 ^ 

g(t) = / /exp j [ (wq + o)d - m)T + yyx 

-T/2 -00 

dtodr 



+ ( ^o - 

2y 



+ w t ] 



T/2 ^ 

= / exp j [cJqT + (OdT + J]iT^ 

-T/2 



+ 



2y 



X 



+ 00 2 

fexp j [- WT + + 03t 



— ^ m] dm dr 

y 
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Call V = 



(gjn + MT - yt) 
/2 y 



V 



2 _ tJJp ' 



2y 



1 2 1 2 
y yx + OJqT + J 



- yxt - Wq! 



Multiplying g(t) by exp -v^ and exp v" 



1 T/ 2 1 

" 2 ^ ^ + 2 1-''^^ - v^] X 

X /^°exp j [w - /2\i v] ^/2 y dudx 

Since ^ ^ 

bince, 2y 2y ^ 2y 

2 . j . Wn 

= - WT + + OJt - — (jJ 

2y y 

If u = - Y - , the second integral will be of the form 

/■2H 

/2]x /^°°exp j du = /Ty /IT exp j 'n/4 



1 2 

which implies g(t) = , exp j (wq! - y yt ■*■ x3 

T/2 



Since 


WqT + wdx + 1 


yx^ - + 


yxt - 


i yt^ + wqP 


which 


implies 


g(t) = 


/■ -T— COS (oJnt " 

\—-) 
y J 


I Wt") 


sin (wd + yt)T/2 




(wd + yt)/2 


which 


implies 


g(t) = 


/,2yT^, , . 

( ^ )cos (Wot - 


I wt") 


sin X 
X 


where 




X = 


(wd + yt)T/2 







X = iihll + ^ 
^22 
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From these results it can be seen that i£ w, = 0, a ^ 

type o£ envelope is obtained with a peak at the origin o£ 

magnitude /^ -ZAcuT ^ . So the higher AtoT the better the output 
ir 

signal. I£ ^ 0 there is a shi£t in the sin x/x curve which 
implies a range error £or doppler £requencies greater than 
zero . 

3 . Matclied-Filter Approach 

The matched- Filter approach is another passive method. 
The analysis o£ the Former will give us a basis For comparison. 
Let the received signal be 

F(t) = cos [(o)ot+ o)(j) t + y yt^] , |tl < T/2^0 elsewhere. 
The output oF the matched Filter will be oF the Form 
g(T) = h (T - 1) dt 

where h(t) =KF(-t) and F[h(t)] = H(u)) is the Filter transFer 
Function. K is a constant Factor to give a unity gain. In 
this case, K = / ^2y^ . So 

g(T,u)d)= ^ itTcos [ (o)o + wd) t + i yt^] 

TT ^ 

cos [u)o(t - t) - iy (t - T)^]dt 



aFter some steps [7] 

g(i^>t^d) = 4"2y\ cos [ (u)^ + ^)'t] sin ^ (T 



cod + yx 



xi)] 



For I X I < T 



or, g(x,Wd) = j (T- |xl)’^cos[(mo + ^^)x] 



Sin X 



, Ixi < T 



where x = (T - |x[) 
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Since x is not a linear function of t, the shape of 
the curve, unlike in the previous case, is not a sin x/x type 
of curve. Also, since does not represent a linear shift 

in the axis due to the non-linearity of t with respect to x, 
if 6Jd 0 the curves become distorted and not symetric v;ith 
respect to the vertical axis as in the previous case (Fig. 7). 
The advantage is of course the maximization of the signal to 
noise ratio. 

4 . Use of Discrete Frequency Sequences in Pulse Compression 
The use of a sequence of frequencies each ^^^ith a pulse 
duration of t in order to get a good autocorrelation function 
is a very common method. 

Various combinations, from linear stepped frequency 
to randomly chosen frequencies, are used in order to increase 
the main lobe with respect to the side lobes in the ambiguity 
function. The basic format is that of Fig. 8 



A 



nr 



- 



f, 


t 










L 


L 



Fig. 8 

where the waveform of the nth segment can be expressed as 
Vj^(t) = A exp j 27r(fj^t + 
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The mathematical analysis for the general case is 
quite difficult and of little interest. Since only a few 
cases are of interest, consider the linear stepped frequency, 
that is, a time function of the form, 

N-1 

v(t) = ^|Q[u(t-nT) - u(t-(n + 1 )t) ] cos + nAm) t 

where = lowest frequency to be transmitted 

Aw = frequency spacing 
N = number of frequencies in transmission. 

Assuming that there is coherency between all frequen- 
cies, that is, cj)^ = 0, vsfhich could be obtained using a frequency 
synthesizer with a master oscillator, the matched filter would 
be of the form indicated in Fig. 9. 




Fig. 9 



39 



ties . 



If xAf > 1, there ivill be significant range ambigui- 

Also, if xAf < 1 the subpulse filters f^ of Fig. 9 

will overlap. That is the reason why in most cases xAf is 

made equal to unity, that is, t = ^- With that assumption, 

the convolution of the impulse response of the matched filter 

h(t) = v(-t) will yield an output 

N-1 

" n=0 ^n ^ (w^ + nAto)t 

where (N-l)x < t < Nx . 

If A =1, for all n, then 
n ’ 

N-1 

e^(t) = exp j w^t ^E^(exp j Awt)^ 

Since I exp j Awt | < 1, then 

e (t) = [exp j [m + (N-1) Am ] t] x sin N (Aa)/2)t 

2 sin (Aw/ 2) t 

Accounting for the autocorrelation of the rectangular 
pulse, the final output time function is 

e (t) = (1 - |xlAf) sin N(Aw/2)t r (N-1) Aw -. 

° sin (Aw/2)t 

Fig. 10 shows the effects of xAf / 1. The nules of 
the output waveform occur at sin NAw/2t = 0. 

Which implies N ^ t = ± m ir 

For N > 50, the first side lobe is 13.46db below the main 
lobe . 
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A very important quantitative factor in all these 
waveforms with sharp autocorrelation function is the pulse 
compression ratio. It is defined as the ratio of peak power 
after compression to peak power before compression. Table II 
gives some data including compression ratios for some typical 
linear pre-pulse generation and compression circuits. 

Another important quantitative factor is the side lobe 
reduction factor which is defined as the power ratio of the 
main lobe with reference to any side lobe. 

Depending on the specific application, it is sometimes 
preferable to sacrifice the main lobe peak power or even width 
(3 db point), but get a higher side lobe reduction. In the 
radar problem those side lobes may be identified with targets, 
the need to suppress them if wanting to use a larger dynamic 
range is obvious. 

There are various types of weighting functions. The 
weighting is applied in the amplitude response of the matched 
filter in order to deliberately produce a mismatch. 

The frequency response of the matched filter is such 
that the signal to noise ratio is maximized. So any mismatch 
will lower that signal to noise ratio. 

The loss in signal to noise ratio due to weighting 
is called the loss factor L^. It is defined as 

T = (S/N) weighte d ^ dtf 

^ (S/N) matched Tf-pW'^ (t) dt 

o)(t) = weighting function 

T = processing time interval or the length of the 
transmitted waveform. 
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Table II 



In the case of discrete frequency amplitude weighting, the 

loss factor is 
N 

[S An] 2 

•-s ' -hs [8) 

N[E am 
1 

The problem of weighting is in effect a way of producing a 
pulse compressed spectrum that yields the wanted waveform. 

In the continuous linear pM case, the spectrum of the 
match filter is of course a pulse in the frequency domain. 
Since the inverse transform of a pulse is a sin x/x type of 
function, that is Achy for the linear FM waveform a sin x/x 
type of output is obtained. From the Fourier transform pair 
it is known that a wide pulse in one domain corresponds to a 
sharp pulse in the opposite domain. So the amplitude weight- 
ing functions that will lower the side lobes will have a bell 
or tapered shape characteristic in the frequency domain. That 
is why the time output of a matclied filter with a pulse in the 
frequency domain as transfer function, is 

1 w 

e^(t) = / a(w) exp j cot dco 

U Z TT - W 



where a(co) is the weighting function; the most conunonly used 
weighting functions are: 

The cosine function where 

a(co) = cos ^ where a(co) = 1 at the center of the 

spectrum . 



1 W 

So, ©nCt) / COS 

U V /TT -w 



IT CO 

2tv 



exp j (0 1 d (0 = — 2" 



2 w 

~ 4 u 7 ^ 



cos cot, 
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The main lobe will be between nules. The nules will 
occur at cos = o 

which implies wt = (2n + 1) , n = 0, ±1, ±2,,.. 

or t = (2n + 1) TT 
2w 

Considering the lobes between nules the side lobe reduction 
is 23.5 db . 

Another class of tapers are the Hamming functions 
which have the general form 

G (w) = a + (1 - a) cos ( -~) 0 < a < 1 

where the center of the spectrum is assumed to be at w = 0. 

In this case the output v;ill be 

= — sin v-'t 1 -a 

^ 2 . wt 2 w ( 1 - 1 ^ ) 

for a = . 54 the side lobe reduction is 42.8 db which represents 
a much larger dynamic range. Table III gives data on side 
lobe suppression on frequency coded waveforms. 

D. MOVING TARGET INDICATOR (MTI) 

The MTI is basically a processor that separates moving 
targets from clutter. Since the clutter spectrum is a very 
low frequency spectrum, the MTI processor must be a higli pass 
filter. In the early MTI processors the main technical 
problems \jere in the processor itself due to instability of 
some elements. Today the main problems result due to the 
existence of a variety of types of clutter whicli make the 
design difficult. If only the elimination of one type of 
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clutter is required, then the efficiency of the MTI is good, 
but if various types of clutter must be handled at the same 
time by the same processor, the efficiency drops. 

Depending on the parameter used, the MTI processors are 
basically divided into three types; those that use the phase 
information of the returned signal, those that use the ampli- 
tude, and those that use both phase and amplitude. The basic 
configuration of the phase processing or coherent MTI is that 
of Fig. 11. Tlie use of two oscillators, the stalo and coho, 
is for up and down converting the waveforms in order to got 
a perfect phase comparison, since at low frequencies the phase 
accuracy measurement increases. Coherency is also obtained. 

Consider two signal returns from the same moving target 
with a separation in time by T = 1/fj^ wliore f^^ = pulse 
rejection frequency. The signals presented to the summer are 

= E sin (o)dt + 4)^) 

E 2 = F sin [ooj(t + T) + (^ 2 'i 

<j>l = (^2 = = 'i’ 

The output of the canceller will be 

Er = - E 2 = 2Esin(^) cos[u)j(t + ^) + (()] 

R 

where it is seen that the cosine waveform is modulated by a 
sin . So, independently of t, there are some doppler 

frequencies where the power out will be zero. These are the 
frequencies where • That is, f^ = nfj^ (multiples 

of the pulse repetition frequency) . 
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Fig. 11 
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Since the doppler frequency is a function of the radial 
2 Vr 

velocity Vj^, f^ = -y— , there are also blind speeds that 
correspond to the blind frequencies (f^ = 

^blind ^ ) , n is an integer. 



Fig. 12 shows the power response for a single delay line 
MTI processor. As it is seen, the multiples of the PRF are 
at nules but the other frequencies are not processed in the 
same way. That is very far away from the ideal case. 

The clutter spectrum beating with the pulse repetition 
frequency is translated to every multiple of the PRF. So if 
the clutter spectrum is, after translated, of the form indi- 
cated in Fig. 13 (a) the ideal MTI processor should have a 
response as in Fig. 13 (b) . 

As it is known the single delay line processor response 
is far away from that of Fig. 13 (b) . One way to improve the 
MTI processor is by using multiple canceller filters. 

For an n-stage cascaded canceller (Fig. 14) 



^So-, ')2n . an.TTfa., 

(^)n =2 sin (— p^) 



[9] 



-R 



where and are the peak output power and peak input power, 
The problem is that if many cascaded blocks are used, the 



(sin function may cut some low velocities of interest. 



2n 



R 



The transmitted waveform has a spectrum as in Fig. 15. 

The fact that there is more than one line in the spectrum is 
of no importance to the MTI processor since all the lines are 
at blind frequencies. 
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Fig. 12 
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Fig. 15 



52 






The amplitude processor, also called noncoherent processor, 
has a block diagram as in Fig. 16. The advantage o£ this 
processor is due to the fact that stability problems almost 
disappear since there is no coherency. Since most high power 
transmitters are not very stable, this is a real advantage. 

The big disadvantage is that it needs a permanent existence 
of a strong clutter return since without it, it won't work. 

Fig. 17 represents the phasor diagram that explains how the 
noncoherent processor works. Since the clutter is of the 
form, = dc value, in a phasor diagram is represented by a 
stationary vector. 

Let Eg(t) be the return of a moving target. 

Then ^ lEg(t)| cos w^t 

Eg(t + T) = cos [w^(t + T)] 

Let ^d^ ~ ‘I’ “ 

to^(t + T) = 4) + A(J)/2 = (J) - Acf)/2 + oj^T 
which implies w^T = Ac() 

Let the total signal at time t (Fig. 17 b) be 
E\(t) = E-; + 

and the total signal at time (t + T) (Fig. 17 c) be 
E^2(t) = E^(t + T) = e"^ + E^^(t + T). 

From Fig. 17 (d) it can be seen that due to E^, there exists 
an amplitude difference between |E^(t)| and |E 2 (t)|. This 
difference in amplitude is going to be used to identify the 
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Fig. 16 
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moving target. In order to get better results, instead of 
the difference in amplitude, the difference of the square of 
the amplitudes will be used. That is the reason for a square 
law detector in Fig. 16. 

From Fig. 17 (b): 

E| = + E| - 2E^E| cos ^ - A0/2) 

(the equation relates only the absolute values of 
the vectors) 

From Fig . 17 (c) : 

E| = E2 - E| - 2E^E^ cos ((J) + A({)/2) 

(only amplitude relations) 

Then E2 = Ej - E^ = 4E^E^ sin (f) sin ^ 

since (p = w^t + A<{)/2 

and A<t> = w ,T = ^ j 1 

R R 

which implies 

E^ = 4E E sin 2-fd_ sin (wt + A<p/2) 
r c s 

So, the blind speeds are the same as in the coherent 
processor, and the transfer function must have the same shape. 
But if E^ = 0, that implies E^ = 0. That is the reason why 
an amplitude processor must have a clutter return, and the 
stronger the clutter return the better since E^ increases 
linearly with E^ . 

Clutter cancellation can also be made at IF frequencies 
instead of at the video part of the receiver (Fig. 18). 
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Fig. 18 



Let the two 
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E 
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2 



where d) 



IF signals separated by T = be 
^ E sin [2TT(f^p ± f^)t - 4 .^] 

E sin [2 tt(£jp ± £^) (t + T) - 4 ,^] 

„R 

IF o 

c 



Er = Ep - Ez = 2E sin [ir (£jp ± £^) T] cos 

[ 2 ir (£jp ± £^) (t + i) - 4 )^] 

is the output o£ the summer. 

The output o£ the phase detector would be 

E^ = E sin [tt(£jp ± £^)T] cos {27r[£^t+ (fjp ± f^) j] ~ 

where the video envelope is o£ the £orm E sin [it (£jp ± £^) T] . 
Since in an MTI processor i£ £^ = 0 the output must be zero 
that means that ■'^fjpT = nir is a necessary condition. That is, 
^IF “ ^ ^ IE Frequency must be a multiple o£ the 

pulse repetition Frequency. 

The existence oF blind speeds can constitute a problem 
since every velocity is not processed in the same way. One 
oF the most oFten used methods to improve that situation is 
with the use oF staggered PRF systems. There are some quanti- 
tative Factors that describe the eFFiciency oF an MTI processor. 
Since they are the descriptors oF the MTI processor, it is 
worthwhile to state their proper deFinitions, 

Improvement Factor is deFined as I = where r^ is 

the output target to clutter ratio and r^ the input target 
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to clutter ratio. This definition reflects the gain as well 
as the clutter rejection of the processor. 

Subclutter visibility is defined as the capacity of a 
radar to detect moving targets in a clutter environment. A 
radar with x db of SCV is one that is able to detect a target 
over a clutter that has a signal x db stronger than the moving 
target. The SCV cannot be used as a parameter for comparison 
between radars since the target to clutter ratio is a function 
of the size of the radar resolution cell. 

The two main problems with MTI processors design are: 
first, a correct modeling of the expected type of clutter, 
and second, the correct shaping of the transfer function 
depending on the type of clutter. In [2] there is an extensive 
description of the different types of clutter and modeling 
processes . 

When the platform of the radar is moving, the stationary 
targets will have a doppler return different from zero. This 
is the typical case of the Airborne MTI pi'ocessors. Various 
techniques of clutter looking and automatic tracking must be 
used. The techniques get much more complicated and the 
processor becomes much more dependent for specific application. 

E. DIGITAL SIGNAL PROCESSING, DFT AND FFT 

Due to high speed, low cost, versatility and almost inde- 
pendence of external condition of recent digital computers, 
the processing of radar signals in digital form covers now a 
wide spectrum of the radar signal processing. 
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The two basic methods of approach are represented in 
block diagram in Fig. 19. In the first case (Fig. 19 a) the 
signal is first sampled and quantized in such a way that the 
signal is converted into a sequence of numbers. The numbers 
are then stored and the processor, using arithmatic and logical 
operations, manipulates the numbers according to an algorithm 
Avhich is a function of the type of filter we want. The 
resultant numbers are then dequantized and passed through a 
low pass filter (D/A) . The second method (Fig. 19 b) is 
basically a frequency domain approach of the problem. The 
signal is quantized and a discrete Fourier transform algorithm 
is applied to the quantized signal. The spectrum of the signal 
is then weighted as a function of the type of filtering re- 
quired. Then, the inverse operations take place. Due to 
recent improvements in faster algorithms to implement the DFT 
tills second method is becoming widely used. The theory behind 
the first method is based on the Z transform theory [10] . The 
variable Z represents a delay and is defined through its Laplace 
transform by 



where 1/T is the sampling frequency. Given an analog function 
f(t), the transform of f(t) will be 



The two basic configurations are recursive and nonrecursive 
filters. In the nonrecurs iv'e case the output is of the form 



Z = exp ST 



F(Z) = Jq f(nT)Z 



E = B. [1 + A, Z 
o 1 1 
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weinhting factors 



where the coefficients are obtained from the coefficients 
of a Fourier series expansion of the frequency domain function 
of the desired filter. 

The feedback or recursive filters are normally obtained 
using the analog transfer function of the required filter and 
substituting S by a chosen relation between S and Z. Usually 
the bilinear transformation 

[11] S = f ^ i-y 

is used since with the correct sampling frequency it is very 
accurate and easy to implement. Figure 20 shows the basic 
configurations for the nonrecursive (Fig. 20 a) and recursive 
(Fig. 20 b) filters. 

In many radar applications the discrete Fourier transform 
is being used especially as a bank of filters for obtaining 
doppler outputs. The understanding of what it is and its 
limitations are thus fundamental. 

In order to get a Fourier transform with a digital pro- 
cessor not only the input waveform must be in digital form 
but also the values of the Fourier transform must be repre- 
sented by discrete values. Figure 21 gives us the graphical 
derivation of the discrete Fourier transform pair. The input 
waveform h(t) (Fig. 21 a) must be sampled in order to be 
manipulated by a digital processor. The sampling operation 
is the same as multiplying h(t) by a string of delta functions 
(Fig. 21 b) . Since multiplication in time is the same as 
convolution in the frequency domain, the resulting spectrum 
is that of Fig. 21 c) . Here it is seen that if the samples 
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are not at a frequency higher than the Nyquist rate there 
will be an aliasing effect that will distort the waveform. 

But since in most systems, and especially in any real time 
radar system, the waveform must be processed in a finite time, 
there are for computation only a finite set of samples. This 
is equivalent to multiplying h(t)A^(t) (Fig. 21 c) by a 
rectangular function (Fig. 21 d) . Again, the previous spec- 
trum must be convolved with the spectrum of the rectangle to 
get the result in Fig. 21 (e) . As expected, the larger the 
, the smaller is the distortion in the spectrum. But as 
stated before, the output must be in a digital format. So it 
is necessary to multiply in the frequency domain |H(f) * AQ(f) 

* x(f)|, by a string of delta functions, separated by l/T^. 

That is equivalent to the convolution in tisne of a string of 
impulses separated by T^ with h(t)A^(t) x(t). The final result 
is in Fig. 17 (g) . As can be seen, the output of the DFT 
processor is not a sampled replica of |H(f)|. Aliasing effects 
and finite time of processing are responsible for the amount 
of distortion introduced. In order to get the mathematical 
expression for the DFT, it is only necessary to make a parallel 
derivation to that of the graphical derivation already dis- 
cussed. The time domain expression for the output (Fig. 17 g) 
will be obtained; then the Fourier coefficients of that 
periodic function represent the DFT. 

Let h(t) be the input function. The operation of sampling 
h(t) is equivalent to a multiplication by 
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+ 0O 



A (t) = E 6(t - KT) 
° K=-oo 



■j" oo 

so h(t)A (t) = E h (KT) 6 (t - KT) 
° K=-~ 



Since >^(t) 



- 4- < t < T - T/2 
2 o 



0 otherwise 



is to be multiplied by h(t)A^(t) 



N-1 

h(t)A^(t)x(t) = E h(KT) 6 (t - KT) 

° K=0 

where NT = T 

o 

then we have to sample the Fourier transform of h(t)A^(t)x(t) 
by a string of impulses in the frequency domain 



Aj(f) ■= £:«(£ - i/T^) 

which infers Aj(t) • F‘ [Aj(f)] - ‘ f 



must be convolved with h (t) A^ (t)x (t) to obtain the time domain 
answer . 

hp(t) = [h(t)A^(t)x(t)J A^(t) = 

N-1 +0O 

= [ E h(KT)6(t - KT)] [T )] 

K=0 

■K» N-1 

= T^ E h(KT)6(t - KT - rT ) ] 

Since hp(t) is a periodic function with a period NT, in 
order to avoid time aliasing x(t) is chosen in such a way that 
the end points of x(t) don't coincide with sampling points. 
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I£ they did coincide, it would generate an additive effect 
at the boundary. But NT is still equal to . Also, the 
Fourier transform of a periodic function is given by the 
Fourier coefficients of a series expansion. So 



where 



H,,(nf^) = = E^'^c 5(f-nf ) 

O'- 0 -’ -«> n ^ o-^ 

c = ^ h^(t) exp -j dt, n = 0, 1,. 

^ o -T/2 ^ ^0 



, T -T/2 

So, c = Y~ f ^ Tq 
^ ^o -T/2 



rLcoK=o ^CKT) 6 (t - KT - rT^) exp 



dt 



Since the integration is only over one period T^ 



To-T/2pN-l 27 rnf 

c = / E 6 (t-KT) exp 

^ -T/2 K--0 o 



dt 



c 



n 



N-1 T -T/2 2 Trnt 

E h(KT) / exp-j--^ 6 (t-KT) dt 

K=0 -T/2 o 



. 2frKnT 

E h(KT) exp -3 — rp 
K=0 



o 



Since Tq ~ NT 



N-1 



n 



E h(KT) exp -3 
K=0 



. 27TKn 



N 



, n = 0 , 1 , 2 , 



and the Fourier transform of hj^(t) is 



Hd C-1^5 ’JZ '/h(KT) exp-j 6 (f - 

K- 0 
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To see that also periodic, it is necessary to 

show that c^^ is periodic Vv^ith a period NT. 



Let n = r + N 



since exp -j 



. 2TrK(r + N) 
^ N 




which implies c^ (r + N) = c^(r) 



which implies (j^) 



so there are only N distinct values that can be evaluated. 
Normally the function 



is called the discrete Fourier transform (DFT) which relates 
N samples in the time domain to N samples in the frequency 
domain. The inverse discrete Fourier transform, is, as expected 



So it is seen that with very simple complex multiplications 
and additions there is conversion from one domain to the other. 

The Fast Fourier transform is an algorithm to implement 
the discrete Fourier transform which, by reducing the number 
of multiplications reduces the processor time to perform the 
algorithm. For a better understanding of the FFT, it must be 
reduced to a matrix form. 

Consider the DFT 




N-1 

E h (KT) exp - j Z-rrnK/N ,n = 0, 1... N-1 
K=0 



1 N-1 n 

h (KT) = i E H exp j 

K=0 




N-1 

x(n) = E X (K) exp -j 
K=0 




0, 1. . . N-1 
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if W = exp -j 2 tt/N, then in a compact form 

n V 

x(n) = W’"'" Xo(K) 



where x(n) = 


x(o) 


, XoCK) = 


>^o(o) 




x(N-l) 




x^(N-l) 




'w° -- 


--- W° 1 






W° -- 


... li(N-l)CN-l) 





the examination of the matrix equation shows that there will 
2 

be N possible complex multiplications and N(N-l) complex 
additions. If N = 2^, and the FFT algorithm [12] is used, the 
number of multiplications will be reduced to Ny/2 and the 
number of additions to Ny . Figure 22 gives a comparison 
between the FFT algorithm and the direct computation with 
respect to the number of required multiplications. 

As it was shown, complex notation is greatly used in the 
mathematical structure of signal processing. So the need to 
represent signals in complex notation becomes a necessity even 
in their implementation. It v\rill be seen how a sine wave type 
of function is converted in order to be represented by a 
complex number. 

In Fig. 23, let the input signal be 

S(t) = A sin [(Wjp ± w^)t + (J)^] 

where w-.,, and w, can be thought as the intermediate frequency 
1 r d 

and doppler frequency in a radar receiver. If the local 
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oscillator (I/Q) generates 2 cos and -2 sin w^pt which 

can be obtained w'ith a 90° phase shifter, tlie in pJiase channel 
(I) and the out of pliase channel (Q) outputs are of the form 

I = A cos (w^t + (j)^) 

Q = A sin (w^t + (|>^) 

Since the low pass filter retains only the difference frequency. 
The I and Q signals can then be thought as the real and imagin- 
ary part of the complex signal 

Z = I + j Q 

If tlie I and the Q channels are sampled and converted to 
sequences of numbers Xj. and y^., the resultant sampled pair 
can be considered a complex digital word. 

\ ^ " |Aj(|exp j <|)j^ 

where = A cos K T + {t>^) 

y = A cos (o) j K T + 4)^) 

|Aj^| = A, ^ <J>o 

assuming a sampling frequency f = ^. Since (J)^, varies linearly 
with K, the complex digitized signal can be interjected as a 
vector of amplitude A that shifts at every sampling time 
KT, by lo^KAT. 



72 







:-iV 




F. SYNTHETIC APERTURE RADAR 



Using the pulse compression techniques already mentioned, 
it is possible to achieve very small values of range resolu- 
tion. In many radar applications such as aerial photography 
not only a very small range resolution is necessary but also 
an equivalent small aximuth resulution is desired. The con- 
ventional method of reducing azimuth resolution by very small 
beamwidths will always require big antennas wliich, especially 
in airborne radars, is a problem since space is a limiting 
factor. Also, the fact that azimuth resulution is an increas- 
ing function with distance makes it many orders of magnitude 
higher than the range resolution even for small distances. 

Consider a radar without a pulse compression processor. 
Then the range resolution 6R is 



6R = 



2 

w 



c = speed of propagation 



B = bandwidth 



If the antenna has an aperture D, the beamwidth is 6g = -jy- 



D ’ 



So the azimuth resolution 6^^ is 



^AZ D 



R = distance to the target 



if B =10 MHZ, f = 5000 MHZ, D = 10 yds and R = 100 mi, 
then ^ 

6^.^ = 130 0 yds 

Synthetic aperture radar is a way of processing radar 
signals in order to get an equivalent antenna aperture much 
bigger than the real antenna to reduce SA^. In the synthetic 
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aperture case a srngle antenna is translated along a line 
The received signals are then stored, and processed after 

the radrating element has travelled a fixed distance. The 

synthetic aperture concent cnn hr. .r- j , 

P an be viewed through two different 

aspects. From a doppler viewpoint or from a i’, 

F j-iiL or trom a linear array view- 
point, the basic principle of the doppler effect viewpoint is 

there exists a one to one correspondence between the 

long track coordinates of a 

oi a reflecting object and the instan- 

tSHGous doppler shift, 

consider Pig. 24 in which an airplane flies at an altitude 
h„ with a constant speed v and with a radar antenna that 
illuminates the area A, The radar antenna does not rotate 
but simply moves with the airplane. 

If there is no pulse compression and the width of the 

transmitted pulse t i-f- a ^ ^ 

p S T, It IS known that the slant resolution 

constant is 



,CT 



R 2 - Tb 



w 



so the corresponding ground resolution constant is 
pR = sec ^ ^ sec ip 

Also, from Fig. 24 it is seen tlmt t - o d i 

seen that L - 3R where 3 is the beam- 

width of the transmitting antenna and is the the along 



track resolution p . 

X 



Since 



3 = 



D 



then 



Px ^eff ~ “d" 
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and as pointed out, it is impossible to increase D beyond 
certain limits in order to decrease p . But if a relation 
is found between doppler frequency and the position x of the 
radar, then it is possible to make an azimuth discrimination. 

From Fig. 25 it is seen that the radial velocity with 
respect to the target is 

= V sin 9 



So the doppler shift will be f , = ^ f^ sin 6 where f is the 
transmitted frequency. But if 6 is small, which is the normal 
physical situation since the total illuminated area is small 
and R is big, then sine=;e = x^-x, which makes 

-c - 2 y . 

d AR ^ o ^ 



So it is seen that with a frequency analysis of the received 
signal, it is possible to get an azimuth resolution that is 
only a function of discriminating between the closeness of 
two frequencies. The basic principle for the explanation of 
synthetic aperture radar using the linear array theory is the 
following : 

In a linear array various elements are fed at the same 
time and the received echo is also received at the same time, 
but due to path differences a receiving pattern is generated. 

In the synthetic aperture case there is only one transmitting 
element that moves with constant speed along a line. So a 
long antenna be formed not by physical means but by signal 

processing. After the radiating element has travelled a dis- 
tance corresponding to the illumination time of the real 
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beam, the stored signals after weighting and phase shifting 
resemble the signals received by a big linear array. 

Two cases arc important to differentiate; the unfocused 
synthetic aperture and the focused synthetic aperture. 

In the focused synthetic aperture radar if 0 = 0 (Fig. 26 a) , 
by signal processing (for R or by real situation (R = «>) 

the signals all arrive with the same phase. In a real linear 
array the angle selectivity is provided only during the recep- 
tion of the signals. But in the synthetic case, since only 
one element is radiating and is moving, it is necessary to 
account for the phase shift due to the transmitted and received 
path. So, due to the fact that all the round-trip phase counts 
for the formation of the receiving pattern, as opposed to the 
real case in which only the received path is responsible for 
the receiving pattern, the effective beamwidth is 



*eff 



2 L 



instead of 



eff 



'eff 



Since the length of the synthetic aperture radar is equal 
to the distance corresponding to the time of illumination, 

Leff = ^ , where D is the horizontal aperture of the physical 
antenna and R the distance to the target (Fig. 24), then the 
azimuth resolution is 

*AZ ^eff ^ 2 ^ 2 

The azimuth resolution is independent of X and R in the 
focused case which is of course an important result. This 
means that the smaller the antenna the higher the azimuth 



resolution. 
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In the unfocused case it is important to account for the 
phase difference between the center and the ends of the equiv- 
alent linear array. A quantitative interpretation of the 
unfocused synthetic aperture radar can be explained as follows 
Consider Fig. 26 (a): Due to the fact that i R, there 

is a phase shift between the center and the end. The total 
voltage is the vector sum of small discrete voltages that are 
received along the path (Fig- 26 b) . If the number of 

discrete voltages is large, it is possible to make the approx- 
imation of Fig. 26 (c) and determine the relationship between 
the resultant voltage V and the linear sum of the discrete 
voltages . If the 3 db point is defined as the boundary 
for the drop in voltage, a must be approximately 7 t/ 2 . But, 
when a - ir/2, the radius r^ (Fig. 26 c) is perpendicular to 
the tangent t^^, implying that the phase difference between 
the voltage at the end and at the center must approximate 45°. 
This sets a limit on So v;hen R = R^ + A/8 (round trip 

< > A/4) the conditions are met. Using Fig. 26 (a): 

(Leff/ 2)2 + R" = (R + X/ 8)2 = r 2 + ^ 

^ 

if is very small compared to the other terms L^££ - / RA 

Since = .455 t / A*~ 

3db 2 

the azimuth resolution would be 6 ^^ “ ^ 

1/2 

Since 6.7 is a function of R , this is a solution between 
the focused case and the conventional case. Figure 27 gives 
a relation between the three cases. 
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The signal processing tlieory for synthetic aperture radar 
can be understood through the analysis of the ambiguity func- 
tion of the received signal. If it is possible to separate 
the range and azimuth as the two ambiguity factors from the 
general ambiguity function, then the determination of the 
azimuth ambiguity function and consequently the azimuth reso- 
lution function arc immediate. 

Consider a transmitted signal f(t). After reflection in 
the ground or target with reflectivity p(x,y,z), the received 
signal will be 

S(t) = Iff p(x,y,z) f (t - ^) dxdydt [13] 

where the integration is made over the illuminated region, 
and R is the distance between (x,y,z) and the radar antenna. 
The processing of the signals in fact recovers a signal pro- 
portional to p(x,y,z). Using matched filter techniques, the 
output will be of the form. 

Op I 

e^Ct) = f S(t) dt 

= ffff p(x,y,z) f (t - f* (t - dtdxdydz 

where R' is the distance from the antenna to the point 
(x',y',z') corresponding to p(x',y',z'). Let's find the 
ambiguity function of f(t). By definition it will be 

i|^(x,y,z,x' ,y ' ,z') = f f(t - f* (t - dt 

where * means complex conjugate. 
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So 

eo(t) = ////ipp (x,y,z) dxdydz, 

is a weighted average of p(x,y,z) where the weighting function 
is i(;(x,y,z,x' ,y' ,z') . 



Let 

then 



£(t) = g(t) exp j (D^t 

_ r 2R^ 2R\ . ,2R 2R', ,, 

= / g(t - — ) g* (t - -^) exp j 0)0 (— - dt 



Assume the energy is transmitted in finite amounts of 
time. If during those finite lengths of time the exponential 
does not vary too much, then during a transmission the expon- 
ential part can be considered constant although varying between 
transmissions. So it can be taken out of the integral. The 
result will be 

i}; = E [ / g(t - g* (t dt] exp -j (-^ 

If gg* is constant during a transmission period, it can be 
moved out of the summation and 

, r/- J U.1 V ■ f2R 2R'i 

ilj = [ f gg^dt] Z exp - — ] 

It is obvious that the first factor is the range ambiguity 
function, and the second, the azimuth ambiguity function. 

So, the azimuth ambiguity function is 

Figure 28 indicates a real situation. From that figure 
it can be seen that: 
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X 



vt 






R' = + fx - X ' J 2 : R + Cx - X ' ) ^ 

o 

if X << and (x^ - x') << , which is basically the same 

restriction given to e in Fig. 25. 

If the radar transmits with a PRR of f„ = 1/T, then x = 

K 

nvT are the positions where transmission occurs. 



which infers 



E exp - j 



x2R 



2R’ 



= ^ exp - 2 ^^ (2xx' -(x*)^) 
o 



= exp + 



a3o(x')^ 



N/2 

E 

-N/2 



exp 



-j 



4tt (x ' /AR ) n vT 

V o 



where L, the 



which infers 



synthetic aperture length is 

L = NvT , V = aircraft speed. 

= pvn i ^ sin (N + 1) 4 ttx ' vT/2ARq 

P ^ R c sin 4ttx ' VT/2XR 

o ' o 



Tx ' 1 ^ 

where exp j - 4 is only a phase factor. The 
o 

sin (N+l)a/sin a is the amplitude factor of interest. The 
azimuth resolution is given by the 3 db points separation. 
If =: 1, the 3 db point is at 



2ttx’L 



and 




1 . 4 X Ro 
ttL 
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Since L = 



6 



AZ 



3R = K R 
o Do 



1.4 

TT 



D 



D 

2 



which is a result in agreement with that already derived. 

In generating the synthetic aperture radar, in order to 
reduce the side lobes of sometimes useful to weight 

the returned signals before combining them. In the case of 
the focused antenna, phase compensation must also be used. 

So, if N returns are processed in order to get 

L = N V T 

the focused processing will be of the form 

Z S [ exp j \p ] W 
n '■ ^ n-* n 

and the unfocused processing of the form 

S S W 
n n 

where is the weighting function and S^ are the discrete 
returns. This is the basic processing technique. Various 
methods, from optical, electronic, as well as acoustic, are 
used to generate the synthetic aperture radar. 
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V . RECENT DEVELOPMENTS IN TWO MAJOR AREAS 



A. DIGITAL MTI 

The structure o£ an MTI processor clearly points into the 
direction of a digital implementation. The increase in digital 
data rates associated with the cost decay of present digital 
components led to the implementation of reliable, flexible, 
low cost digital MTI processors. 

Fig. 29 is a typical block diagram of one of the types of 
MTI digital implementation. I and Q channels are used in 
order not to lose 3 db on the average due to blind phases. 

The structure is parallel to the analog MTI processor. The 
use of A/D's and D/A's, storage devices instead of delay lines, 
make up the differences with respect to the analog processor. 
The capacity of the memory depends basically on the number of 
range cells and on the use of multiple cancellers. 

The quantization process introduces a new dimension to the 
noise problem: that is, the quantization noise. Quantization 
noise is present along all the dynamic range of the processor, 
but the errors at the extremes are of particular importance 
to the MTI processor. At the lower level, clutter cannot be 
cancelled below one bit of quantization j at the upper level, 
abrupt clipping distorts the signal creating additional noise. 
This extra noise affects the improvement factor (I) of the 
processor. If it is assumed that signal and quantization 
noise are uncorrelated, then for a single canceller, the 
improvement factor is [14], 
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where 




Pg, correlation coefficient of quantization error for 
clutter alone 

correlation coefficient of quantization error when 
both clutter and signal are present 



- noise power 
E = signal voltage 

= input power clutter 

P(T) = correlation coefficient of clutter signals. 



Since a ^ /E^ 0 and p (T) = 0 

e e 



1 - P(T) + OeVPi, 



which differs from the equivalent analog factor only by the 
term a VPi^- ^ 

Also, for a double canceller, the improvement factor 
is [14] 

1 



I. = 



1 - i p(T) H- I p(2T) + a^VPic 



For the triple and higher order cancellers, the expressions 
get more complicated but the derivation is the same for all 
cases. Also in I^, the only difference from the analog 
expression is the factor Assuming a uniform error 

distribution 

P 2 

a " = 

^ 12 ( 2 ^'- 1 / 2) 2 

where E,^ = saturation voltage of the digital register 
n = number of bits used in the quantization. 
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is called the improvement factor limitation. 

Figure 30 relates the improvement factor limitation with the 

total number of bits for different values of P. /E 2 

1 C m 

In many applications the fact that the MTI processor has 
blind speeds at multiples of the PRF is completely undesirable. 
The method usually used to attenuate this effect makes use of 
staggered PRF. The scan-to-scan staggered PRF is less efficient 
but requires less complex hardware than the pulse- to- pulse 
stagger, but for some applications the first one is inadequate. 
The staggered PRF technique is based on the use of different 
spacing betv/een transmitted pulses which, at the receiver, are 
properly delayed (de- staggered) , such that the spacings at 
the input of the delay line cancellers are all equal. The 
necessary use of delays makes staggered PRF perfectly matclied 
for digital techniques. Since the concept is completely 
defined, the recent trend has been focused in the development 
of optimization techniques in order to improve the performance 
of the processor. The first difficulty with any optimization 
problem is, of course, the definition of optimum. Some try 
to optimize the improvement factor, others to optimize tlie 
signal to clutter gain (SCG) within a finite number of fre- 
quency slots, or even the optimization of some defined indicator. 

An important factor in the optimization process is the 
statistical model tliat is chosen for the clutter power dis- 
tribution. It has been verified that a Gaussian distribution 
centered at zero doppler frequency is adequate in most cases 
as a power density function of clutter. Various standard 
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I 



deviations have been calculated for different types of 
clutter [ 2] . 

Figure 31 is a block diagram of a staggered PRF processor. 
The pulse staggered sequence is such that the i^^^ pulse is 
delayed (i - 1)T + AT^ seconds with respect to tlie first pulse. 
In the processor, after de- staggering , the pulses are weighted 
(W^ ' s in Fig. 31). By observation, the impulse response of 
the processor is 

h(t) = 6(t - AT^) + IV 2 6[t - (T + AT 2 )] 

+ ... + Wj^ 6{t - [(N-l)T + ATj^]} 
which implies a transfer function 

H(oa) = W. exp {-jo) [(i - 1)T + AT.]} 

and a filter power response 

G(to) = H(o)) H*(w) , * E complex conjugate. 

With G(w) and the power spectral density of the clutter, 
different optimization criteria can now be devised. One factor 
is common to all techniques. The A^'s, W^'s and the clutter 
variance, assuming a Gaussian distribution, are the parameters 
that will optimize the processor, given an optimization pro- 
cedure. In the analysis, it is always assumed a step scanned 
array antenna in order that for each position, N pulses are 
transmitted. Figure 31 is a processor for a given range cell. 
If there are K large cells, the general block diagram is 
presented on Fig. 32, v^^here ' s and the summing network is 
presented for the K. cell. 
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The ideal MTI filter is a high-pass filter with rejection 
from zero frequency to the maximum frequency of the clutter. 
Rejection may also occur after the maximum desired doppler 
frequency but it is not absolutely necessary. 

For a better understanding of staggered PRF some general 
concepts should be explored. If a constant PRF is used, not 
only the dc component of the doppler spectrum but the whole 
clutter spectrum is translated to multiples of the PRF. It 
is shown [15] that with staggered PRF the clutter spectrum 
still is translated to multiples of the average PRF, but the 
spectrum of the signal is dispersed into N (number of A^^’s) 
separate frequency lines. So, if a target produces a doppler 
frequency that is a multiple of the basic PRF, only one part 
in the whole signal is distorted by the clutter spectrum 
where the remaining M-1 parts of the signal may possess enough 
power to be detected. An MTI filter which notches in a region 
around dc as well as in multiples of PRF can eliminate the 
clutter, but if staggered PRF is used only one part in N of 
the signal is lost, thus increasing the probability of detec- 
tion of targets with potential blind velocities. 

In [16] an optimization process, based on a performance 
index P defined as 

^u 1 

P = /. u df 

^1 [S2(f)]^‘ 

is developed. and are the lower and upper bounds of the 

velocity region, S^(f) is the signal to clutter gain and M is 
a parameter that reflects the emphasis the criteria puts on 
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spectral amplitude variations. As M increases, smaller 
spectral amplitudes have more influence in the final result 
of the integral. 

For an understanding of the significance of P, it should 
be noted that good clutter attenuation as well as a high im- 
provement factor correspond to a good performance index. 

Using the optimization procedure described in [16] and using 
M = 8, practical results (figures 32, 33, 34) for different 
values of N and (B^ - B^) were obtained. It can be seen that 
as N increases, the minimum value for the signal- to-clutter 
gain in the velocity region increases as well as the improve- 
ment factor I = SCR(f) (Fig. 32, 33), but with a fixed N and 
the velocity region increased on the upper bound, the minimum 
value of SCR(f) in the velocity region decreases (Fig. 33, 34). 

As was said earlier, the clutter variance a is also an 
important parameter in the optimization process. Fig. 35 
relates a with the improvement factor for the case of N = 5 
and B^ = 10 B^. In the solid line the interpulsc periods 
T + were optimized for the value of a = 0.02 B^ and the 

coefficients were optimized in parallel with the variation 
of a in the horizontal axis. On the dotted line, both coef- 
ficients and interpulses were optimized for a = 0.02 B^, and 
the variations of I with a were recorded. This situation 
occurs when the real clutter parameters are not those used 
in the design of the filter. It is apparent that there are 
not too many differences between the two cases. 
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SCR gain es a function of Doppler frequency; 
^ = = 0.02/?, . Optimum intervals: 1 .1 00. 

1.094, 1.000, Optimum coefficients: 1.000 —3 003 
3.155,-1.152. 



Fig, 32 




SCR gain as a function of Doppler frequency; 
A/ = 5. D^, = 1 05j . a = 0.02£l, . Optimum intervals: 1.111, 
1.000, 1.091, 1.058. Optimum coefficients: 1.000, 

-4.258, 6.242, -4.031, 1 .048. 



Fig. 33 



SCR gain as a function of Doppler frequency; 
N ^ — 40^, . o —0.07By . Optimum intervals: 1 .007, 

1.021, 1.000, 1.299. Optimum coefficients: 1.000, 

—3.784, 5.432, -3.273, 0.626. 




Fig. 34 



Improvement factor as a function of clutter spectral width 
o; // = 5, = 1 05, . Coefficients for dotted curve (optimun^ for a = 

0.025,): 1.000, -^.258, 6.242, -4.031, 1.048. Intervals for both 
K)lid and dotted curves (optimum for o ~ 0.025,): 1.111, 1.000, 
1.091, 1.058 {c.f.. Fig. 8). 




Fig. 35 
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Another criterion to optimize the staggered PRF processor 
is developed in [17]. Instead of maximizing the average sig- 
nal- to-clutter gain over the entire velocity region, the 
s ignal- to - clutter gain is maximized in small Af intervals 
v\fithin the velocity region. The number of intervals is, of 
course, directly related to the complexity of the processor. 

Using the optimization algorithm developed in [19] , a final 
block diagram for the processor is reached (Fig. 36, 37). In 
Fig. 36 the input signal after being split in I and Q channels 
is weighted with the optimum set of weights A and B for each 
channel. The results are then combined in order to get an out- 
put of positive doppler as well as an output of negative doppler. 
Fig. 37 represents the implementation of the filter that gener- 
ates the sets of A and B weights for each of the M frequency 
sub - intervals in a channel. The T^'s are the interpulse delays 
and the 

, 0<i<N-l 

a . . b . . 

0<j<M-l 

are the coefficients determined by the optimization process. 

It should be noted that in this optimization process the 
interpulse periods are not optimized: only the Aveights are 
optimized. Fig. 38 relates the signal- to-clutter gain aver- 
aged over the entire velocity region with the signal-to-gain 
maximized in each of eight frequency sub- intervals . Both 
N = 5 and N = 10 are presented. A substantial improvement 
due to optimization is evident. Another optimization tech- 
nique is developed in [18] . 
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Fig. 37 
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The optimized processors are called Constrained Improvement 
MTI Processors (CIP) . Given a specified improvement factor, 
the number of pulses to be processed and the PRF stagger 
sequence, the mean square deviation from a constant response 
in the velocity region is minimized. Two main approaches can 
be given to the problem. In the first, the PRF stagger 
sequence is fixed and the weighting coefficients are chosen 
in order to optimize the processor; in the second, the weight- 
ing coefficients are fixed while the stagger sequence is chosen 
in order to optimize the processor. Defining 

f = ^ , F' = maximum PRF 

Fr 

F^ = unstaggered PRF 

figures 59, 40 and 41 represent the optimum filter response 
with a four pulse return, f = 8, a = 0.01 and I = 30 db, 
using respectively a linear PRF stagger with ±20% interpulse 
variation, and a sinosoidal PRF stagger with ±10% (Fig. 40) 
and ±90% (Fig. 41) interpulse variation. As it can be seen, 
as the variation increases the response becomes more uniform. 
Comparing figures 42, 39 and 43, in which the number of pulses 
are respectively 3, 4 and 6 for a common interpulse variation 
of 20%, it can also be seen that the response improves with 
the number of pulses processed. 

A comparison between optimization procedures, different 
from a direct comparison between responses, can be devised 
using the fraction of frequencies for which a response is 
less than some specified value as a comparison parameter. 

If the 0 db scale corresponds to the mean value of the 
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Frequency response of four-pulse CIP using sinusoidal PRF 
stagger with ±10 percent interpulse period variation. 

S 




Fig. 40 




Freguency re$pon:e of four-pulse CIP using sinusoidal PRF 
stagger with t 90 percent interpulse period variation. 



Fig. 41 
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Frequency response of three-pulse CIP with ± 20 percent 
interpulse period variation. 



Fig. 42 



Response of six-pulse CIP using sinusoidal PRF staggar 
with ± 20 percent variation in interpulse period. 







Fig. 43 
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frequencies of interest, that means that 50% of the frequen- 
cies are above as well as below the 0 db level. This leads 
to a cumulative distribution function in which instead of a 
fraction of frequencies, the term probability is used. Fig. 

44 is a cumulative distribution of the responses of a four 
pulse CIP and a four pulse processor using the Ref. [16] 
optimization criteria. This comparison criteria gives a 
better performance for the CIP processor since its curve is 
closer to being the ideal step response at 0 db than the other 
one, but on the average they are similar. 

^ Another type of digital MTI processor is called the matrix 
MTI . In this case there is no parallel between the analog and 
the matrix MTI processor. This processor (Fig. 45), in addi- 
tion to the tremendous flexibility in the modification of the 
shape of the response curve, has associated with it a very 
simple clutter locking mechanism that shifts the response 
curve wlienever there is an average velocity associated with 
clutter. To simplify the explanation of the processor in Fig. 
45, only two levels of quantization are used. Basically the 
processor is a coherent system that compares the phase differ- 
ence between two successive returns and weights that difference 
in phase according to the shape it is intended to give to the 
response curve. Since dj) = w^dt, by v\feighting d4> in effect 
the doppler is being weighted. 

The phase angle of a signal return is defined by 

4) = the ^ Y 
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Fig 




and the amplitude by (I^ + Q^) .J Since in the case of 
Fig. 45 there is only one quantization level for 1 and Q, 
that means that the only possible phases for I are 0° and 
180° and for Q, 90° and 270°, assuming a signal to noise ratio 
much greater than 1. When the average phase difference is 0°, 
it implies a stationary target, but if d4>,p = 180° that indicates 
an optimum velocity target. Since each vector return can be 
located in all of the four positions ±I - jQ, there are 4x4 
possible combinations for two consecutive returns. So, what 
the compression matrix block in Fig. 45 does is to generate 
signals proportiona], to the phase difference output d4>.j. 
between two consecutive vector returns, whicli in this partic- 
ular case can take the values 0°, 90°, 180° and 270°. Each 
of the four outputs has a corresponding weighting factor W 
and tlie result is summed to generate the MTl output. Fig. 46 
indicates the four types with the vector combinations that 
generated them, as well as the weighting factors used in order 
to approximate an jMTI response curve. The final response curve 
corresponds to a statistical average of when uniformly 

random phases are introduced at the input. In order to have 
a signal proportional to the clutter velocity, it is only 
necessary to subtract the outputs of the region (d) (Fig. 46) 
from those of region (c) and divide by the number of range 
gates, that is, for two quantization levels and a uniformly 
distributed phase input 

dd = rcd) - (c)U9o°) 

^ N 

and VAR d<i> = — (45°) 

3N 
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Fig. 46 
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So, in a large clutter environment if N is large, can 

smoothly adjust the local oscillator to place the MTI null 
at the mean value of the clutter doppler. 

Consider now an n quantization bit circuit. For each 
channel there are 2^ possible quantization levels. If only 
one quadrant of possible phase differences is considered, it 
can be seen that there are (Fig. 47) 

2 ^ 2 ^ 2 ^ ^ , 

— X — - — + 1 

2 2 2 



possible phases. The -2^ ^ + 1 parcel results from the fact 
that of all the diagonal combinations (45° zone) only one 
combination can be counted. As an example, let's use two 
quantization bits. That means that there are 



4[2 



2\\-2 



2H-1^ 



1 ] 



= 12 
n = 2 



possible phases in all four quadrants, which also implies 12 
possible phase differences. 

Figure 48 (a) represents the response and figure 48 (b) 
is the matrix output (phase) for each 12 x 12 possible com- 
bination of two consecutive vector returns. Since there is 
flexibility in modifying the weighting coefficients, it is 
possible to adapt the filter as a function of the type of 
expected clutter spectrum. Again, the mean clutter velocity 
can be determined by subtracting the outputs of the (b) region 
(Fig, 48) from those of the (m) regions and applying the 
result to the VCO associated with the local oscillator. 
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Total number of combinations per quadrant 
counting the diagonal combinations only as 
one = 



^n-1 «n-l 

2 X 2 



+ 1 



Total number of combinations for the four 



quadrants = 



4 [2 



n-1 



>n- 1 



1 ] 



= 2^" - 2"*1 + 4 



Fig. 47 
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B. DIGITAL SAR 



The synthetic aperture radar (SAR) is a typical example 
of signal processing where digital techniques are being sub- 
stituted for the former optical processors. With the present 
technology mini- computers can compensate for the aircraft 
movement, and, as with a real antenna, the synthetic antenna 
can be steered and even scanned. The basic problem with SAR 
processors is the need for high storage capabilities and very 
fast data rates. The film as an optical storage device became 
almost impractical since a real time display is impossible. 

The use of storage tubes is also inadequate due to low effi- 
ciency and poor dynamic range and stability. Witli the present 
high speed A/D converters (100 MHZ), low cost, and high speed 
compact digital storage devices, it becomes feasible for real 
time SAR processors. 

Table IV shows a time overview of the late achievements 
in SAR techniques. 

Figure 49 is a block diagram of a digital synthetic 
aperture radar. Theoretically, the processor must, for each 
range gate, perform the integral 

/ s(t)r(t)dt 

AT 

where s(t) is the signal return and r(t) the correlator 
reference function. 

r(t) = Aj^(t) exp -j 4)j^(t) 

Aj^(t) = weighting function to control the side lobes 
of the synthetic antenna pattern. 

(j>R(t) = phase reference that tracks the phase of s(t). 
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Date 


" Development 


1951 


Carl-Wiley postulates doppler beam 
sharpening concept. 


1952 


University of Illinois demonstrates 
beam sharpening concept. 


1957 


First SAR imagery using optical 
correlator is produced. 


Mid 1960 's 


Analog electronic SAR correlation 
demonstrated in non real time. 


Late 1960's 


Digital electronic SAR correlation 
demonstrated in non real time. 


Early 1970's 


Realtime digital SAR demonstrated 
with motion compensation. 



Table IV 
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Since (j)j^(t) follows the phase of s(t) but with an unknown 
difference, the signals are processed in an I (in phase) and 
Q (90*^ out of phase) channel in order to prevent signal losses. 
The digital processor operates in the following way: The radar 
returns pass an A/D converter that samples the signal at a 
rate, at least higher than the Nyquist rate, and separates 
the returns in range bins which are approximately equal to 
the range resolution. The number of bits used per range bin 
are a function of the desired dynamic range of the processor. 
The digital data passes through a buffer and prefilter which 
translates the A/D rate to the rate of the correlator. The 
bulk memory stores the data corresponding to an integration 
time, AT. The data is then correlated with the reference 
generator output in order to form an azimuth line with a 
length corresponding to the number of range gates. After 
correlation, new data enters the memory and the process re- 
peats itself to produce a nev.’’ azimuth line. 

From Fig. 50 it can be seen that the bulk memory in bits 
must be equal to the number of cells times the average number 
of bits per cell, so 

BM = bulk memory = 2K N N„ 

Si K 

ar 

Nj^ = number of range gates = — 

^r 

A 'T’ 

N = number of azimuth data lines = 

The factor of 2 is because of the use of an 1 and Q channel. 

But if, for example, the desired range resolution is .30 m, 
that corresponds to a bandwidth 



112 



synthetic antenna 
pattern 







113 



I 





BW = _ 3 X 10^ m/sec 

2 X 30 m 



5.0 MHZ , 



Avhich implies an A/D sampling rate per channel of 5.0 MHZ. 

On the other hand, for a 1 KHZ PRF the unambiguous range is 
80 mi, which for a typical AR = 10 mi of mapping, means that 
there is still an equivalent 70 mi excess time to process that 
data. This results that in fact there is no need for the 
correlator to work at such high data rates. Using a buffer 
that accepts the data at a very high data rate during a short 
period of time but delivers it at a lower rate to the corre- 
lator, it is possible to reduce the correlator rates by many 
orders of magnitude. In fact, analyzing the function of the 
correlator, it can be seen that it must in AT seconds produce 
a number of outputs equal to the total number of cells, that 
is, (L/R^)Nj^, where L/r^ = CA is the azimuth compression rate. 
Since the real antenna illuminates (L/r )Np cells and is ready 
to receive new data after T seconds, the correlator rate (CR) 
is 





T 



defining 



Np = CA = total number of doppler filters 



K = azimuth over sample factor ^ 1 



e 



Then 




if Np = 500 and f =1 KHZ for one doppler channel CR = 



.5 MHZ < A/D sampling rate. 
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On the other hand, if the spacing corresponding to V/T 



is much finer than the range resolution r^, that means excessive 
data is being used. If the PRF cannot be lowered due to power 
or doppler considerations, then a prefilter can be used to 
effectively reduce the PRf = f , by a factor f /f where f 
is the prefilter sampling rate. The prefilter will reduce 
the correlator rate by a factor f^/f^. If the required reduc- 
tion is in order to have V/T of the order of r^, then 
V L 1 ^ F 

f ci — =: - — Tp - — rp , which implies a new correlator rate, 

s r r r T ’ ^ 

a a 

CR = K K N2 f 
os s F R r 

K = prefilter over sample factor 
os ^ ^ 



= synthetic array weighting 



These two factors, the bulk memory and the correlator rate, 
determine the type of design for a SAR processor. The main 
objective is to reduce both CR and BM. Today's technology is 
attacking this problem in two ways: devising correlator algor- 
ithms which will reduce the bulk memory and arithmetic, at low 
cost, with fewer power consuming elements. The first approach 
is being made by parallel and series combinations of correlator 
channels and prefilters, as well as with the use of FFT algor- 
ithms. Another method is the use of pulse compression tech- 
niques but processed only after the SAR processor. This will 
reduce which will reduce BM. The hardware improvements 
are in CCD memories and LSI at much lower costs. 

As it was said before, with today's digital techniques 
various mapping modes are possible as well as motion compensation. 
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Figure 51 shows the four types of mapping used today. Tlie 
Squint mode (Fig. 51 B) is only a variation in angle of the 
side-looking SAR. The doppler beam sharpened mapping is a 
PPI representation with an equivalent antenna of very high 
azimuth resolution. The spotlight is a mapping Avhere a liigh 
resolution snapshot map is generated. Motion Compensation 
includes various functions: clutter tracking, focusing accel- 
eration compensation and real antenna stabilization. Fig. 52 
is a block diagram of a general SAR processor. The blocks 
are basically the same as those of Fig. 49 except for the 
motion compensation blocks which are now included. 

To achieve lower data rates, various types of algorithms 
are presently used. Lower rates will be achieved at the 
expense of more complex hardware. The general concept tliat 
applies to all of them is the reduction of rates by the in- 
crease in the number of doppler filters. 

Figure 53 is a block diagram of the multi-channel prefilter 
processor approach [19] where m is the number of channels. 

The prefilter rate becomes 



and since the number of filters Np is reduced byl/m, the total 
arithmetic rate (TAR) is 



R. 



P 



m Nj, f 
R r 



TAR = R + R = aini+^ 
c p ^ m 



a = f 
1 R r 



a = K K Np Np/AT 
2 OS s F R' 




a 



m 



a 



0 
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tA) SIDE LOOK STRIP MAPPING 




(C) DOPPLER BEAM SHARPENED 
MAPPING 




(B) SQUINT K50DE STRIP MAPPING 



A 




A 



(D) SI'OTLICHT MAPPING 



Fig. 51 




Fig. 52 
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implies that the value o£ m that minimizes TAR is 

m = [K K N^,/£ 

'■ os s h r ■' 

and the minimum TAR is 

TAR (9Lu)'4 = 2(a a 

min 

= 2(K K f /aT)^/2N„i\’ 

The two-stage correlator (Fig. 54) is another algorithm [19]. 
The idea is again frequency domain division but only with one 
prefilter. The equivalent correlator rate will be the sum of 
the two correlator rates. So 

TAR = Kg Nj^ Np/AT) M 

+ K Ri N„/MAT 

os s R F 

where M is the number of first-stage correlation channels. 

1 / 

As in the previous case, there is a value of M = (Np) ^ that 
minimizes TAR. With the FFT algorithm, it is possible to 
reduce even more the total rates. Table V is a summary of 
the calculated values [19] of bulk memory as well as total 
arithmetic rates for the SAR signal processing algorithms. 

Another approach, completely different from those already 
presented to solve the problem of high data rates on the 
correlator (Fig. 49) of the SAR processor, is through the use 
of parallel processing using associative memory [20]. Fig. 55 
represents a block diagram of the associative memory and the 
related input/output registers. In the associative memory. 



119 




o . 

to 



•S o 



c o 
o o 
a o 

^ Q* 
to V. 

c o 

^ i 
o Z 

n 

c *5 

pj 

'H 

g 3 

■n w 
•r: to 
o 2 
^ o 

W v^ 

x; o 

to 2 
C to 
o t> 

E c 

P3 •'^ 



.S:P J3 

r: 

> .ti 

H o 

> to 

"1*^ 
a o 

S *5 

15 

o 2 
*o i::: 

Ho 

UP 



G 

rt M 

^ O 

O g 

o tli 

C .s 

C t) 

S5 

{^.S 



> 

0 

r—^ 

rO 

oj 

H 



120 



ADDRESS 
• I 

REGISTER 






INPUT 



1 




SEARCH /V 





L 



MASK REGISTER 

zzzTcr: 



ASSOCIATIVE 

MEMORY 



READ REGISTER 



TAG 
REGISTER 



OUTPUT 



Fig. 55 



121 



multi addressing, multi read/write, multi logical and arith- 
metic operations are performed. The search/write and the 
mask registers arc used to generate proper codes for multi 
accessing. The operations are done in a bit serial basis but 
on all words in parallel at the same time. Ellis [20] showed 
that for an X band radar with range and azimuth resolutions 
of 10 meters, 1000 range cells, maximum range of 150 km and 
maximum aircraft velocity of 200 m/sec, a 50 msec correlation 
time is necessary: this implies a 40 nsec multiplication time. 
Using associative memory, it is possible [20] to perform tlie 
operations of a correlation period in 31.48 msec. The advan- 
tages of liardware parallel operations are thus obvious over 
the conventional process when speed is an important factor. 
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VI. CONCLUSIONS. TRENDS 



It was shown that due to improvements in integrated 
circuitry, better microprocessors, faster DFT/FFT algorithms 
and lower cost of digital logic, most analog processors are 
being replaced by digital processors. Very high data rates 
continue to be the major problem of digital processors, but 
with today's techniques of parallel processing, general pur- 
pose digital radar signal processors are already in use. 
Despite the quantization noise inherent to digital processors, 
good reliability, extreme flexibility and low cost of digital 
processors still give them a tremendous advantage over the 
former analog processors. 

The future demand for more reliable and sophisticated 
radar systems will be a function of cost and military neces- 
sities. If the cost of digital logic trends lower as is its 
present trend, and military demands continues high, digital 
processors will play an increasingly important role in modern 
radar system's implementation. 
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